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Abstract. Based on the daily data of American and Chinese stock markets, the 
dynamic behavior of a financial network with static and dynamic thresholds is 
investigated. Compared with the static threshold, the dynamic threshold suppresses 
the large fluctuation induced by the cross-correlation of individual stock prices, and 
leads to a stable topological structure in the dynamic evolution. Long-range time- 
correlations are revealed for the average clustering coefficient, average degree and cross- 
correlation of degrees. The dynamic network shows a two-peak behavior in the degree 
distribution. 
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1. Introduction 

A financial market is a complex system composed of many interacting units, which 
exhibits various collective behaviors pfl |2j [3J 111 El El [7]. How to extract its structure 
information has attracted much attention of physicists, and remains challenging. For 
example, the hierarchical structure of financial markets has been investigated with the 
minimal spanning tree and its variants [H El QUI [Til H2] • With the random matrix theory, 
business sectors and unusual sectors can be clarified, and topology communities are also 
revealed HJ El [151 EU EH EB] • 

The complexity theory provides a powerful tool to understand complex networks. 
In recent years much progress has been achieved on the fields from biology to sociology 
[T9l |20| [2T| 122] . Different models and theoretical approaches have been developed to 
investigate the formation of network topologies, the functions emerging from networks, 
the combined mechanism of topologies and functions, etc [231 1211 [251 EH1 127] . Most 
efforts of these schemes focus on the collective properties of the system in a steady state. 
However, complex systems such as the financial markets are essentially nonstationary 
in time. For each time step, the system shows a particular topology, induced by the 
cross-correlations of individual stock prices. The topology dynamics is important for the 
full understanding of the network structure. However, it is rarely touched or elaborated 
in detail, to the best of our knowledge. 

So far, static topological properties of financial markets from the view of complex 
networks have been widely investigated [281 ESI EDI EH EH [33]. An ordinary way to 
construct a financial network is to take individual stocks as nodes, and set a threshold 
to create edges. If the cross-correlation between two stocks is larger (smaller) than 
the threshold, then link (cut) an edge. A number of contributions have adopted an 
artificially given threshold for static cross-correlations. However, it does not extract the 
dynamic characteristics of the stock markets [29 1 [32 1 133]. For example, an outstanding 
company with a great number of correlated companies is driven bankrupt by an extreme 
event, and all its linkages are then cut down. This may lead to a temporal variation 
of the topological structure of the market. Therefore, a real financial network should 
be dynamic, to capture the dynamic evolution of the topological structure. On the 
other hand, as is well known, the volatilities of the stock prices exhibit a fat-tailed 
probability distribution. Large fluctuations of the volatilities may greatly influence the 
network structure and network stability, and therefore should be taken into account in 
understanding the topology dynamics. 

In this article, we examine a dynamic financial network based on the American and 
Chinese stock markets. For a comparative study, both static and dynamic thresholds are 
respectively adopted in the network construction. Our purpose is to study the statistical 
properties of the dynamic network, and more importantly, to investigate the temporal 
correlations of the topology time series, such as the time series of the average clustering 
coefficient, the average degree and the cross- correlation of degrees, by applying the 
detrended fluctuation analysis (DFA). Special attention is put on the dynamic effect of 
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the thresholds on the network structure and network stability. 

In section II, we present the data analyzed and construct the dynamic network with 
static and dynamic thresholds. In section III, we investigate the time-correlations of the 
average clustering coecient, the average degree and the cross-correlation of degrees, and 
examine the degree distribution of the network. Finally, the conclusion comes in section 
IV. 



2. Dynamic financial network 

For a comprehensive understanding, we analyze two different stock markets, the New 
York Stock Exchange (NYSE) and the Chinese Stock Market (CSM), representing the 
mature and the emerging markets, respectively. For both markets, we investigate the 
daily data of 259 individual stocks, with 2981 data points from the year 1997 to 2008 
for the NYSE, and 2633 data points from the year 1997 to 2007 for the CSM. 
We define the price return 

Ri(t,At) = ]nPi(t + At) -\nPi(t), (1) 

where P%{t) is the price of stock i at time t, and the time interval is set to At = 1 day 
in this article. For comparison of different stocks, we normalize the price return to 

n(t) = ^ !— , (2) 



where o~i = J < R\ > — < Ri > 2 is the standard deviation of Ri, and < . . . > is the 
time average over t. We define an instantaneous equal-time cross-correlation between 
two stocks by 

Gij(t) = ri(t) rj (t), (3) 

Our financial network is constructed in the following way: take individual stocks as 
nodes and set a threshold ( to create edges. At each time step, if the cross-correlation 
Gij{t) > £, then add an edge between stocks i and j\ otherwise, cut the edge. Due 
to the dynamic evolution of Gij(t), the connections between stocks vary with time, 
regardless of the static or the dynamic threshold. This results in a dynamic topology 
of the financial network. To understand the robustness and stability of the network 
structure, introducing a proper threshold is very important. We first consider a static 
threshold ( ~ Q s , 



N N T 



ft ^S£l c » (l1 (4) 

where T is the total time interval, and iV is the total stock number. Since Q s is 
the average cross-correlation over the total time and all stocks, the static threshold 
is reasonable in the sense of providing a uniform standard for cross-correlations of all 
stocks. As shown in figure HJa) and (c), Q s takes rather small values, 0.17 for the NYSE 
and 0.37 for the CSM. However, the cross-correlation Gij(t) is defined by the individual 
price returns, and fluctuate according to the price dynamics. The static threshold may 
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suffer from the large fluctuation of the cross-correlation Gij(t). Hence, we introduce a 
dynamic threshold ( ~ Qd(t), 

9 N N 

QM = N(N _ n E E gm, (5) 

V / 1=1 j=i+l 

<5d(t) takes only the average over all stocks in a single time step, and therefore fluctuates 
synchronously with the cross-correlations Gij(t). As shown in figure Q](b) and (d), large 
values of Qaif) are observed for both the NYSE and the CSM, in addition to the small 
values for most time steps. Therefore, the dynamic threshold may suppress the large 
fluctuation induced by Gij(t), and accordingly creates a stable network structure. In 
this paper, we consider the static threshold from ( = 0.25Q S to 4Q S , and the dynamic 
threshold from ( = 0.25Qd to 6Qd- 

In figure [2J the average clustering coefficient C(t) is displayed for ( = Q s and 
C = Qd(t). The average clustering coefficient C(t) is defined as 



where is the clustering coefficient of node % [23], denoting the ratio of the triangle- 
connection number of node % to the maximum possible triangle-connection number of the 
node. It is observed in figure [2] that the average clustering coefficient C{t) of £ = Qd(t) 
is apparently more stable than that of ( = Q s for both the NYSE and the CSM. The 
common characteristic shared by £ = Q s and ( = Qd{t) is that their average clustering 
coefficients are large. Especially for the dynamic threshold ( = Qd(t), the average 
clustering coefficient C(t) keeps its values around 0.88 and 0.85 respectively for the 
NYSE and the CSM at a highly clustering level, implying a close relationship of stocks 
in the financial markets. 

3. Topology dynamics 

Topological properties of a complex network are usually described by the average 
clustering coefficient, the average degree and the cross-correlation of degrees. 
We investigate the topology dynamics by computing the time-correlations. The 
autocorrelation function is widely adopted to measure the time-correlation. However, 
it shows large fluctuations for nonstationary time series. Therefore, we apply the DFA 
method (SUES]. 

For a time series A(t'), we eliminate the average trend from the time series by 



introducing B{t') = YJ [A(t") — A ave ], where A ave is the average of A(t') in the total 
t"=i 

time interval [1,T]. Uniformly dividing [1,T] into windows of size t, and fitting B(t') to 
a linear function B t (t') in each window, we define the DFA function as 



iV 8=1 



(6) 



= k k E W) - B t(t')} 2 , 



(7) 



t'=l 
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In general, F(t) will obey a power-law scaling behavior F(t) ~ t e . The exponent > 1.0, 
0.5 < 9 < 1.0 and < 9 < 0.5 indicate an unstable, long-range correlating and anti- 
correlating time series, respectively. 9 = 0.5 corresponds to the Gaussian white noise, 
while 9 = 1.0 represents the 1/f noise. 

3.1. Average clustering coefficient 

The DFA function of the average clustering coefficient C(t) is shown in figure EJ For 
the static threshold ( = Q s , a two-stage scaling behavior is observed with a crossover 
phenomenon in between. For the NYSE, the exponent 9 takes the values 9 = 0.76 
for t < t c and 9 = 1.04 for t > t c . For the CSM, 9 = 0.71 for t < t c and 9 = 0.97 for 
t > t c . The crossover time t c ~ 25 days. This result indicates that the average clustering 
coefficients are temporally correlated for t < t c , then transit to the 1/f noise for t > t c . 
In contrast to that of ( = Q s , the DFA function of ( = Qd{t) shows a clean power-law 
behavior, with the exponent 9 = 0.59 for the NYSE and 0.60 for the CSM, indicating 
the long-range time correlation. 

To investigate the robustness of the above results, and the stability of the network 
structure during the dynamic evolution, we may adjust the level of the thresholds. For 
the static threshold, we consider two alternatives ( = 0.5Q S and ( = 2Q S . As shown 
in figure EJ^a) and (c), one does not find any power-law behavior. In other words, the 
network structure is rather sensitive to the specific value of the static threshold. For 
the dynamic threshold, we observe that the DFA function of C (t) remains qualitatively 
the same for ( > Qdif). In figure [3] (b) and (d), for example, the results of ( = 2Q d and 
3Qd are displayed. A clean power-law behavior is detected with the exponent 9 = 0.61 
and 0.61 for the NYSE, and 9 = 0.62 and 0.64 for the CSM. Higher dynamic thresholds 
yield stronger time correlations (i.e., larger values of 9). This implies that the edges 
created by large Gij(t) are rather stable. 

To further understand the topological pattern, we calculate the time- averaging 
clustering coefficient C for different thresholds, 



As shown in the inner panel of figure E](a) and (c), C of the static threshold ( < 0.5Q S 
is close to 1, but as the threshold increases to 1.5Q S , C sharply drops to 0. As shown in 
the inner panel of figure [3]^b) and (d), the time-averaging clustering coefficient C shows 
a mild decay for the dynamic threshold, and presents still rather high clustering even 
for the threshold ( = GQ d (t) for both the NYSE and the CSM. 

3.2. Average degree 

The average degree K(t) of the network is defined as 



1 t=\ 



(8) 



K{t) = -Y,Ht), 

i=l 



(9) 
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where fcj(t) is the degree of node i, denoting the number of the nodes directly connected 
with it. For £ = Q s , Kit) exhibits large fluctuations between and N — 1, as shown 
in figure H](a) and (c). K(t) = N — 1 indicates that every node directly connects to all 
other nodes in the network, while K(t) = corresponds a set of isolated nodes. The 
fluctuation of the CSM is obviously stronger, and K{t) often comes rather close to and 
N—l. For ( = Qd(t), Kit) is almost bounded by two envelope curves, as shown in figure 
@Jb) and (d). The upper envelope is around Kit) ~ 130 to 140, indicating that one node 
is directly connected to about one half of the nodes on average, while the lower envelope 
is around Kit) ~ 70 to 80, indicating that one node is directly connected to about 
one quarter of the nodes on average. The upper envelope implies a high-low symmetric 
distribution of the cross-correlations Gij{t), with one half of the data points below the 
average and the other half of the data points above the average. The upper envelope 
sounds reasonable, for the larger Gij(t) are usually well above the average. However, 
the lower envelope suggests a high-low asymmetric distribution of the cross-correlations 
Gij(t). Why there is such a lower envelope remains to be understood. 

We then compute the DFA function of the average degree Kit). Similar as that of 
the average clustering coefficient, unstable behavior is detected for the static threshold, 
whereas stable network structure is observed for the dynamic threshold, as shown in 
figure [51 For the static threshold, the DFA function for the NYSE shows a power-law 
behavior with the exponent 9 = 0.61 for ( = 0.5Q S , a two-stage scaling for ( = Q s , and 
no power- law behavior for ( = 2Q S . The DFA function for the CSM shows a power-law 
behavior with the exponent 9 = 0.65 and 0.71 for ( = 0.5Q S and Q s , and no power-law 
behavior for ( = 2Q S . These results can also be understood from the time- averaging 
degree 

K=^EK(t), (10) 
1 t=i 

As shown in the inner panel of figure 0(a) and (c), the network has a large number of 
edges for lower static thresholds, however, the number of edges sharply falls off to nearly 
zero as the threshold increases up to 1.5Q S . 

For the dynamic threshold, a robust power-law behavior is observed. The exponent 
9 of ( = Qd{t) is measured to be 0.58 for the NYSE and 0.56 for the CSM respectively, 
somewhat close to 9 = 1/2 of the Gaussian white noises. But as ( increase, the network 
structure stabilizes, with the exponent 9 = 0.60 for the NYSE and 0.61 for the CSM. 
The time-averaging degree K is also found to decay slowly as the dynamic threshold 
increases, and the network shows significant connections for a high threshold as QQd{t), 
as shown in the inner panel of figure [5](b) and (d). 

Why is the dynamic threshold crucial in the analysis of the network structure of 
financial markets? One important reason is that the volatilities fluctuate strongly in 
the dynamic evolution. It induces large temporal fluctuations of the cross-correlations 
of price returns, as is characterized by Qd{t) in figure [TJ Thus the static threshold 
leads to dramatic changes in the topological structure of the network. However, the 
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dynamic threshold proportional to Qa(t) suppresses such a kind of fluctuations, and 
results in a stable topological structure of the network. For example, we have calculated 
the time-averaging degree in two typical periods of time for the CSM, i.e., t = 1550 to 
1650 with small volatilities and t = 2530 to 2630 with large volatilities. For the static 
threshold ( = Q s , the time- averaging degree is 39 and 98 respectively, far away from 
the time- averaging degree K = 68 in the total time interval. For the dynamic threshold 
C = Qd(t), however, it is 107 and 104 respectively, both around K = 107. Similar results 
are obtained for the NYSE, and also for the clustering coefficient C(t). 

3.3. Degree distribution 

The degree distribution function P(k) describes the heterogeneous properties of nodes. 
To obtain a better statistics, we take the degrees of all time steps as an ensemble. Figure 
[6] shows the degree distribution of the static and dynamic thresholds for the NYSE and 
the CSM. For both thresholds, it is observed that the degree distribution is different 
not only from the Poisson distribution, but also from the power-law distribution of a 
scale-free network. For the static threshold, the degree distribution is rather sensitive to 
the threshold value, and the form of the degree distribution changes dramatically. For 
the dynamic threshold, a two-peak behavior is observed with one peak around k = 
and the other peak around k = 140 ~ 150 for both the NYSE and CSM. The peak at 
k = simply tells that there are a large number of isolated nodes in the network. The 
peak at k — 140 ~ 150 indicates that one node is directly connected to about one half of 
the nodes in the network, i.e., the cross-correlation distribution is hight-low symmetric 
for the dynamic threshold. 

The two-peak structure of the degree distribution for the dynamic threshold 
explains the upper envelope of the average degree K(t) in figure HI Moreover, the 
features of the degree distribution is rather robust for different dynamic threshold values. 
It indicates that the large cross-correlations are much bigger than the average cross- 
correlation, so that certain increase of ( does not alter the network topology. 

3.4- Cross-correlation of degrees 

The so-called assortative or disassortative mixing on the degrees refers to the cross- 
correlation of degrees [351 EZ] • The " assortative mixing" means that high-degree nodes 
tend to directly connect with high-degree nodes, while the "disassortative mixing" 
indicates that high-degree nodes prefer to directly connect with low-degree nodes. The 
cross-correlation of degrees is defined as 

( ) M' 1 E a j a k a - [M~ l J^a \(Ja + k a )f 

n) M-^ a \Ul + ki)-[M-^ a \U« + k a )r [ } 

where j a , k a are the degrees of the nodes at both ends of the a t h edge, with a=l,...,M. 
At a certain time, r > 0, r = and r < represent the assortative mixing, no assortative 
mixing and disassortative mixing, respectively. 
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In figure r(t) of ( = Q s and Q d {t) is shown for the NYSE and the CSM. r{t) 
fluctuates between the interval [—1,1], flipping between the assortative mixing and 
disassortative mixing during the time evolution. Defining the time-averaging cross- 
correlation of the degrees, 

r=^2>(f), (12) 
1 t=i 

we obtain r = 0.00 and -0.20 of the static threshold for the NYSE and the CSM, 
r = 0.36 and 0.22 of the dynamic threshold for the NYSE and the CSM. In other 
words, the cross- correlation of the static threshold shows no assortative mixing or 
the disassortative mixing, while that of the dynamic threshold exhibits the assortative 
mixing. 

To study the memory effect, we again compute the DFA function of r(t). As shown 
in figure the DFA function of £ = Q s shows a power-law behavior with the exponent 
9 = 0.73 for the NYSE and 0.77 for the CSM. For ( = 0.5Q S , no power-law behavior 
is observed. It further confirms the unstable network structure for the static threshold. 
For the dynamic threshold, the DFA function of r(t) shows nearly a same power-law 
behavior for ( = Qd(t), 2Qd(t) and 3Qd(t), with the exponent 6 = 0.60 for the NYSE 
and 0.63 for the CSM. 

The community or sector structure identifies different interconnected subsets of 
networks [13j [HJ EE EEJ [39]. To further understand the topological structure of the 
financial networks, we may investigate the dynamic effect of economic sectors. The 
random matrix theory is a representative approach to such a problem, e.g., one may 
analyze the eigenvalues and eigenvectors of the cross-correlation matrix of price returns 
[T3l HU H5] . With a similar procedure, we first introduce the normalized individual 
degrees ki(t) = < k { >)/cr fc ., with being the standard deviation of fcj. We then 
construct the cross-correlation matrix F of individual degrees ki(t), whose elements 

F lj = ^j:h(t)~k j (t), (13) 

1 t=l 

and compute its eigenvalues and eigenvectors. For the network of the NYSE, we do 
observe that the largest eigenvalue Ao corresponds to the market mode, while other 
large eigenvalues A, typically represent different economic sectors, as shown in figure El 
For the dynamic threshold, the results are robust when ( changes from Qd{t) to 6Qd{t). 
For the static threshold, however, reasonable results are obtained only around ( = Q s . 
For the CSM, the dynamic effect of the standard economic sectors is weak, and we need 
a careful analysis to reveal its unusual sectors as reported in Refs. [T3] , 

4. Conclusion 

We investigate the topology dynamics of a financial network by a comparative study with 
static and dynamic thresholds, based on the daily data of the American and Chinese 
stock markets. For both stock markets, the dynamic threshold properly suppresses 
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the large fluctuation induced by the cross-correlations of individual stock prices, and 
creates a rather robust and stable network structure during the dynamic evolution, in 
comparison to the static threshold. Long-range time-correlations are revealed for the 
average clustering coefficient, the average degree and the cross-correlation of degrees. 

The average clustering coefficient and average degree for both static and dynamic 
thresholds are large, indicating the strong interactions between stocks in financial 
markets. A two-peak behavior is observed in the degree distribution for the dynamic 
threshold, very different from the power-law behavior of a scale-free network. 
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Figure 1. Q s and Q d (t) for the NYSE and CSM are shown in (a), (b), (c) and (d). 
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Figure 2. The average clustering coefficients C(t) corresponding to Q s and Qd{t) for 
the NYSE and CSM are displayed in (a), (b), (c) and (d). The time-averaging values 
of C{t) are 0.78, 0.88, 0.68 and 0.85 respectively. 
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Figure 3. The DFA function of the average clustering coefficient C(i) is plotted on 
a log-log scale. The inner panel shows the time-averaging clustering coefficient for 
different threshold values. Dashed lines are the power-law fits, (a) and (c) are of the 
static thresholds for the NYSE and the CSM. Circles, crosses and triangles are for 
C = 0.5Q S , Q s and 2Q S , respectively, (b) and (d) are of the dynamic thresholds for 
the NYSE and the CSM. Circles, crosses and triangles are for £ = Qd{t), 2Qd(t) and 
3Qd(t), respectively. 
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Figure 4. The average degree K{t) is displayed, (a) and (c) are of the static threshold 
C = Q s for the NYSE and the CSM. (b) and (d) arc of the dynamic threshold ( = Q d (t) 
for the NYSE and the CSM. 
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Figure 5. The DFA function of the average degree K{t) is plotted on a log- log scale. 
Inner panels show the time-averaging degrees for different threshold values. Dashed 
lines are the power-law fits, (a) and (c) are of the static thresholds for the NYSE and 
the CSM. Circles, crosses and triangles are for ( = 0.5Q S , Q s and 2Q S , respectively, 
(b) and (d) are of the dynamic thresholds for the NYSE and the CSM. Circles, crosses 
and triangles are for £ = Qd(t), 2Qd{t) and 3Qd(t), respectively. 
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Figure 6. The degree distribution P(k) is displayed on a log-log scale, (a) and (c) 
are of the static thresholds for the NYSE and the CSM. Circles and crosses are for 
C = 0.5Q S and Q s respectively, (b) and (d) are of the dynamic thresholds for the 
NYSE and the CSM. Circles, crosses and triangles are for £ = Qd{t), 2Qd(t) and 
3Qd(t), respectively. 




Figure 7. The cross-correlation r(t) of degrees is displayed, (a) and (c) are of the 
static threshold ( = Q s for the NYSE and the CSM. (b) and (d) are of the dynamic 
threshold ( — Qd(t) for the NYSE and the CSM. The time-averaging values of r(t) are 
0.00, 0.36, -0.20 and 0.22 respectively in (a), (b), (c) and (d). 
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Figure 8. The DFA function of the cross-correlation of degrees is plotted on a log-log 
scale. Dashed lines are the power-law fits, (a) and (c) are of the static thresholds for 
the NYSE and the CSM. Circles and triangles are for £ = 0.5<5s and Q s respectively, 
(b) and (d) are of the dynamic thresholds for the NYSE and the CSM. Circles, crosses 
and triangles are for £ = Qd(t), 2Qd(t) and 3Qd(t), respectively. 
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Figure 9. The absolute values of the components Ui of stock i for the first four largest 
eigenvalues of the cross-correlation matrix F are displayed for the NYSE. Stocks are 
arranged according to economic sectors separated by dashed lines. A: Basic Materials; 
B: Conglomerates; C: Consumer Goods; D: Finance; E: Healthcare; F: Industrial 
Goods; G: Services; H: Technology; I: Utilities, (a) is of the static threshold £ = Q s , 
and (b) is of the dynamic threshold £ = Qd(t). 



